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Abstract 

Detailed comparisons are made between observations 
of scattered light upstream of the head of the HH 1 jet 
and predictions of simple scattering models. It is shown 
that, in order to unambiguously determine the velocity 
of the head of the jet (bow shock) with respect to the 
upstream dust, existing spectroscopic observations are in- 
sufficient and that spectropolarimetric observations of the 
scattered light are necessary. Such an independent mea- 
sure of the bow shock velocity is important in order to 
test "multiple outflow" theories of Herbig-Haro jets. It 
is also shown that the scattering dust must have a very 
forward-throwing scattering phase function ((cos 9) ~ 0.7) 
and slight evidence is found for a dust-gas ratio that is 
higher than average. 

Subject headings: ISM: dust, extinction — ISM: jets and 
outflows — ISM: reflection nebulae — polarization — stars: 
individual (HH 1) 

1 Introduction 

Herbig-Haro (HH) objects are apparently the result of the 
interaction of high velocity outflows from young stars (or 
from protostellar objects) with the surrounding environ- 
ment (see, e. g., the review of Reipurth 1992, which has a 
very extensive list of references). It has recently been re- 
alized that some of these objects can be quite successfully 
modelled as bowshocks, which could be formed, e. g., at 
the head of a jet-like flow ejected by the outflow source 
(the first detailed HH object bowshock models were dis- 
cussed by Hartmann k Raymond 1984). 

Even though at least partial agreement between bow- 
shock models and observations has been found for many 
HH objects (see, e. g., Hartigan, Raymond k Hartmann 
1987), detailed comparisons between predictions and ob- 
servations including line profiles, line ratios, and narrow 
band images have apparently only been carried out for the 
high excitation objects HH 1 (Raga et al. 1988; Noriega- 
Crespo, Bohm k Raga 1990; Noriega-Crespo, Bohm k 
Raga 1989) and HH 34 (Morse et al. 1992), and the low ex- 
citation object HH 7 (Curiel 1992). For HH 1 and HH 34, 
the best agreement with the observed spatially resolved 
line profiles and line ratios is found for bowshock models 
with velocities in the i>bs ~ 150-200 km s _1 range (Raga 
et al. 1988; Morse et al. 1992). On the other hand, the 
low excitation object HH 7 is better reproduced with a 
t; B s ~ 100 km s" 1 model (Smith 1991; Curiel 1992) 

Another observational constraint is provided by the fact 



that the observed proper motions and radial velocities al- 
low us to directly determine the true spatial velocity of 
HH 1 (v s « 352 km s- 1 ; Herbig k Jones 1981) and HH 34 
(v s k, 320 km s" 1 ; Heathcote k Reipurth 1992; Eisloffel 
k Mundt 1992). In both cases, the spatial velocities of 
the sources are likely to be at least an order of magni- 
tude smaller (e. g., no proper motion has been detected 
for the source of HH 34, and its radial velocity is a factor 
~ 5 lower than the one of the HH object), so that the 
vs values also approximately correspond to the velocities 
of the HH objects relative to the outflow sources. Unfor- 
tunately, no proper motion has been detected for HH 7 
(Herbig k Jones 1983), so that its spatial velocity is not 
known. Because of this, in the following discussion we will 
concentrate on the objects HH 1 and HH 34. 

From the results described above, it is evident that for 
HH 1 and HH 34 a clear discrepancy exists between the 
i>bs ~ 150-200 km s _1 velocities deduced from the bow- 
shock models and the vs ~ 300-350 km s _1 spatial veloci- 
ties. There are three plausible explanations for this quite 
remarkable discrepancy (see, e. g., the discussion of Raga 
k Kofman 1992) : 

1. Due to the difficulty of measuring the positions of the 
diffuse HH knots, the true errors in the proper motion 
measurements are much larger than the estimated er- 
rors. 

2. There is some fundamental problem in the bowshock 
models (e. g., the geometry assumed for the bowshock 
is incorrect). 

3. The environment of the HH objects is moving radially 
away from the source at a velocity i>e ~ 150kms _1 . 

For the case of HH 1, new determinations of the proper 
motions (Raga, Barnes k Mateo 1990) clearly confirm the 
previously determined values (Herbig k Jones 1981), so 
that, at least for this object, the first explanation appears 
to be unlikely. Also, as we have described above, the in- 
tensity maps and the spatially resolved line profiles and 
line ratios of HH 1 and HH 34 are well reproduced by 
bowshock models. Because of this, it would be somewhat 
surprising if the i>bs /vs discrepancy were the result of a 
fundamental problem with the bowshock models. 

We are then left with the third explanation. In this 
scenario, the bowshocks are moving into a medium that 
is expanding radially away from the source with a veloc- 
ity i>e ~ 150kms _1 . In this way, a bowshock with a 
spatial velocity vs ~ 350kms _1 has a relative velocity 
i>bs ~ 200kms _1 with respect to the preshock gas, in 
good agreement with the measurements. However, this 
explanation has the somewhat undesirable characteristic 
of introducing a large velocity for the environment as an 
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extra free parameter of the model. If such an environmen- 
tal velocity indeed existed, it would either imply that the 
(optically detected) outflow is moving into a slower (un- 
detected and possibly less well collimated) outflow (see, e. 
g., Stocke et al. 1988), or that it is moving into the slower 
"tail" of a previous "outflow episode" which was ejected in 
the past by the same source (see, e. g., Raga et al. 1990). 

This conclusion appears to be favored by a new observa- 
tional result. Solf & Bohm (1991) (hereafter SB) observed 
faint, highly blue-shifted line profiles in the near environ- 
ment farther away from the source than HH 1. This is 
in sharp contrast to the very low radial velocities mea- 
sured in HH 1 itself (a direct result of the fact that the 
HH 1/2 outflow axis is almost on the plane of the sky). 
Noriega-Crespo, Calvet k B6hm (1991) (hereafter NCB) 
have interpreted these highly blue-shifted line profiles as 
the result of scattering of the emission of HH 1 by dust 
present in the near environment of the object. These au- 
thors find that in order to reproduce the observed line 
profiles, it is necessary to assume a relative velocity be- 
tween HH 1 and its environment of k, 200 km s -1 . This 
measurement would be completely consistent with the sce- 
nario described above (in point number 3). 

However, this apparent confirmation of the existence of 
a moving environment in HH 1 depends on the identifi- 
cation of the blue-shifted line profiles (observed close to 
HH 1, see SB) as HH 1 emission scattered on environmen- 
tal dust. Polarization measurements would clearly confirm 
if this is indeed the case. 

In Henney (1994) (Paper I), a general theoretical frame- 
work was constructed for the interpretation of spectropo- 
larimetric observations of dust-scattered emission lines. 
Predictions were presented of the spatially- and velocity- 
resolved intensity and polarization that would result from 
the Rayleigh scattering of the emission of a moving object 
from surrounding dust. The details of the extension of 
these models to non-Rayleigh scattering from larger dust 
grains will be presented in Henney, Axon & Raga (1994) 
(Paper III). In this paper, such models are fitted in de- 
tail to observations of the upstream light in the HH 1 jet 
in an attempt to extract information about the physical 
parameters of the jet and its surroundings. 

2 Scattering Models for HH 1 

2.1 Integrated Scattered Light 

Solf & Bohm (1991) (SB) present a contour image of the 
HH 1/2 complex in the light of [Sn]AA 6716, 6730 (their 
Figure 1). From this image, it is possible to very crudely 
estimate the scattering optical depth of the surrounding 
dust cloud in the following fashion. If it is assumed that 



the lowest two contours in the region around and upstream 
of the leading condensation HH IF are entirely due to 
scattered light and that all the higher contours are due to 
intrinsic emission, then, by measuring the area under each 
contour, it can be calculated that the ratio of scattered to 
intrinsic flux is ~ 0.2. For a spherical cloud, as is shown in 
Paper I, this ratio is equal to the scattering optical depth 
to the source, so 

T scat = a H n H (R c -Ro)~0.l7, (1) 

where «h is the hydrogen number density, cr H is the mean 
scattering cross-section per hydrogen atom at the appro- 
priate wavelength and R c , Ro are the radii of the scat- 
tering cloud and the source. The value of <r H depends on 
the dust-gas ratio and on the composition and size distri- 
bution of the dust. Possible variations in this parameter 
are discussed in Section 2.2 below but initially a value of 
4.6xl0 _22 cm 2 is adopted, as by NCB, to best allow com- 
parison with their results. For the same reason, a value 
for the distance to HH 1 of 500 pc, and for the radius of 
the source (working surface) of 3" are used, giving 

(i?i5-0.2)n 75 o~4.4, (2) 

where is the scattering cloud radius in units of 15" 
and «75o is the cloud's hydrogen number density in units 
of 750 cm -3 . The radius and density are expressed in these 
units because they are the values derived by NCB in fit- 
ting their radiative transfer model to the longslit spec- 
trum of SB. Their alternative model has R c = 28.5" and 
n H = 500 cm -3 , so they find the right-hand-side of the 
above equation to be 0.8 or 1.13. The sole reason for the 
discrepancy between these values and that of the current 
work is that NCB estimate a much smaller optical depth 
(r scat = 0.03-0.04) and it is therefore worth investigating 
the validity of the assumptions behind the derivation of 
equation (1). 

Inaccuracies can arise either because of errors in deter- 
mining the ratio of scattered to intrinsic flux or because 
the geometry and physical state of the scattering cloud 
render equation (1) inappropriate. Considering the flux 
ratio first, it is unlikely that a significant proportion of the 
light in the lowest two contours of SB's image is due to 
intrinsic emission. This is because the [Sll] longslit spec- 
trogram along p. a. 329° shows that the faint upstream 
emission is kinematically distinct from both the emission 
from the working surface of the jet and that from the fore- 
ground nebula (unlike in the case of Ha, where there is 
significant nebular contamination). Much more likely is 
that a significant amount of the scattered light falls in 
the region of the higher contours and is thus swamped 
by the intrinsic emission. This is expected on theoreti- 
cal grounds since the models suggest that the intensity of 
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scattered light should be strongly peaked near the source, 
especially for forward-throwing scattering (see Paper III), 
and evidence for it can be seen in the broad wings of the 
line profile from the working surface. Hence, the amount 
of scattered light is probably underestimated rather than 
overestimated, so equation (1) is a lower limit on r scat . 

Considering next the idealizations behind equation (1), 
it is possible that extinction should be included since, with 
a mean grain albedo u> ~ 0.5 at this wavelength (Draine 
& Lee 1984), the extinction optical depth would be ~ 0.4. 
The apparent intensity of the source I' is then given by 



I' = I exp {-r ext } = I exp {-r s 



(3) 



since any light absorbed or scattered will no longer appear 
to come from the source. The scattered light, however, 
will only be attenuated by absorption since any further 
scattering out of the beam will be exactly compensated by 
scattering into it. On the other hand, the path travelled in 
the cloud by a photon which has been scattered once will 
be longer than that of an unscattered photon by a factor 



(3 = R{1- 



-{l-R 



2 sin 2 1 
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where 9 is the scattering angle and R is the radius at which 
the photon was scattered, in units of R c . The mean value 
of (3 for the scattered light 1 can then be calculated as 
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(3X(6)R 2 sin6» d6 dR, 



(5) 



where X(6) is the scattering phase function. For small 
dust grains, this will be the Rayleigh phase function (see 
Paper I), and, for larger grains, a modification (West 1991) 
of the empirical Henyey-Greenstein (H-G) phase function 
(Henyey & Greenstein 1941). This latter function has one 
free parameter g, which describes the asymmetry of the 
scattering and is related to the mean cosine of the scat- 
tered light via 

~3(4+£ 2 ) 



(cos6»} = (/- 



(6) 



5(2+ff 2 )' 

It is preferable to use an empirical phase function, rather 
than one derived from simple physical theory, when con- 
sidering the angular dependence of scattering from large 
dust grains because the scattering properties are signifi- 
cantly affected by the porosity and surface roughness of 
the grains, which cannot easily be modelled in a straight- 
forward fashion. This topic will be discussed further in 
Paper III of this series. The quantity (/?) is plotted in 



1 Since this is only a first order correction to the scattered inten- 
sity, it is not necessary to consider multiply-scattered photons 
when calculating this path. 
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Figure 1: Differential extinction of source and scattered light 
is a result of two competing factors: a singly scattered pho- 
ton will, on average, have travelled a distance in the dust 
cloud longer than that of an unscattered photon by a factor 
(/?) (solid line) but, on the other hand, only absorption need 
be considered in calculating the attenuation of scattered light 
(since out- and in-scattering will balance) whereas the source 
light is attenuated by absorption and scattering. The result- 
ing differential extinction factor (dashed and dotted lines) is 
shown for a range of albedos u> as a function of the asymme- 
try parameter of the scattering phase function. 



Figure 1 as a function of scattering asymmetry parame- 
ter g (solid line) for the modified H-G phase function (cf. 
West 1991 and Paper III). The apparent intensity of the 
scattered light I' s is then related to the scattered intensity 
if there were no extinction I s by 



I' S = I S exp {-(l-w)(/3)r scat /w} 



(7) 



and so the ratio of apparent scattered intensity to source 
intensity is 

■77 = ~r exp {-[(1 - w)(/3)-l]r scat /w} = T scat exp {(5r scat } , 
-'o -'o 



where a relative extinction factor 



(8) 



(9) 



has been introduced. This relative extinction factor is also 
plotted as a function of scattering asymmetry parameter in 
Figure 1, for various values of the mean grain albedo u>. As 
can be seen, if the albedo and the scattering asymmetry 
are both small, then 6 is negligible and the source and 
scattered light are both attenuated to the same degree, so 
the extinction can be safely ignored. Even with relatively 
high albedos and a very forward-throwing phase function, 
the extinction factor is always less than unity, as can be 
seen by considering the limit of equation (9) as (/?) tends to 
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Figure 2: Boost factor e for scattered intensity from a prolate 
ellipsoidal scattering cloud over that which would be seen 
from a spherical cloud having the same optical depth to the 
source. The major axis of the cloud lies in the plane of 
the sky and e is plotted against the asymmetry factor g of 
the scattering phase function (modified H-G) for a variety of 
cloud eccentricities. 



one. In this worst case of 6=1, then, using equation (8), a 
new estimate for the scattering optical depth of r scat ~ 0.15 
can be made. Notice that this is not significantly smaller 
than the original estimate. 

Turning now to the geometry of the scattering cloud, 
equation (1) can be extended to non-spherical clouds by 
replacing r scat (the optical depth along the line of sight to 
the source) with a mean optical depth, averaged over all 
scattered photons, which is given by 



error by less than a factor of two. It must be remembered, 
however, that R c in that equation is the cloud radius along 
the line of sight to the source, which is the minor axis 
of the ellipse in the example considered above. The ma- 
jor axis of the ellipse (which will determine the maximum 
extent of scattered light upstream of the source) could, 
depending on the eccentricity, be much greater. A non- 
spherical cloud would also affect the value of (/?) but this 
will not be calculated here. 

Combining the above considerations and adopting g = 
0.6, ui = 0.55 and e = 0.9 (giving 6 = 0.85, e = 1.1) leads 
to a new lower bound on the scattering optical depth to 
the source of r scat > 0.13. While this is lower than the 
original estimate, it is still much higher than the value 
obtained by NCB. Obviously, there are further ways in 
which the analysis could be refined. For instance, none 
of the above calculations take account of the finite size of 
the source, despite the fact that this is implicit in equa- 
tion (2). However, the crudity of the determination of 
the scattered light fraction means that further refinement 
would be unwarranted. 

The extinction optical depth to HH 1 can be deter- 
mined independently, from the observations of Solf, Bohm 
& Raga (1988), who measure a reddening E[B — V] = 0.43 
from which (assuming a total-to-selective extinction ratio 
R v = 3.1; Mathis 1990) it follows that r ext = 1.22. If all 
the reddening is local to HH 1, then (with u> = 0.55) this 
implies r scat = 0.67 but, since the reddening could occur 
anywhere along the line of sight, this is an upper limit to 
the scattering optical depth and is hence quite consistent 
with the value derived above. 
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cll.o.s. 



is the cloud radius along the line of sight 
to the source and e is the ratio of scattered intensity to 
that which would be seen from a spherical cloud having 
the same scattering optical depth to the source. For the 
simple case of a prolate ellipsoid of eccentricity e whose 
major axis lies in the plane of the sky, then e is given by 
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This is plotted in Figure 2 as a function of scattering 
asymmetry parameter for eccentricities of 0.3-0.99 and 
it can be seen that, even for highly eccentric ellipsoids 
(e = 0.99 corresponds to a major:minor axis ratio a/&~7), 
the enhancement to the scattered intensity is not large, 
especially for high values of g. Hence, it seems that, even 
with a very non-spherical cloud, equation (1) will be in 



Spatial Distribution of the Scattered 
Light 



The [Sll] spectrogram of SB is taken with the slit lying 
along the HH 1 jet and so corresponds to Aperture A de- 
scribed in Paper I. It is possible to estimate the velocity- 
integrated scattered intensity as a function of distance 
along this slit. The slope of the scattered brightness pro- 
file measured in this manner is an independent observa- 
tion from that of the total scattered intensity and should, 
in principle, allow the radius of the scattering cloud to 
be determined. Unfortunately, the situation is rather less 
straightforward than this and the derived cloud radius is 
strongly dependent on the assumed cloud shape and scat- 
tering phase function. In Figure 3, three different models 
are presented, all of which fit the brightness profile equally 
well. 

Model 1 is for a spherical cloud with Rayleigh scattering. 
The surface brightness profile is plotted for scattering from 
a finite source (solid line) and a point source (dot-dashed 
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Spherical Cloud, R ( = 18", g = 0.0 



Spherical Cloud, R ( =44.5", g =0.6 



TABLE 1 

Scattering Models that Fit the Observed 
Brightness Profile 
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Figure 3: Three models for the brightness profile in Aper- 
ture A. Measurements from the [Sll] spectrogram of SB at 
2" intervals are shown as filled circles. The edge of the source 
is marked by a vertical line and the source profile (smeared by 
seeing of 1") by the dotted line. a. Model 1. The solid line 
shows the profile of the model (source plus scattered light) 
and the dot-dashed line shows the result of calculating the 
scattering using a point source. The dashed and dot-dot- 
dot-dashed lines show the degree of polarization for these 
two cases, b. Model 2 (non-Rayleigh scattering). The solid 
line shows the profile of the model (source plus scattered 
light), which uses a scattering asymmetry parameter § = 0.6, 
and the dot-dashed line shows the profile of Model 1 (§ = 0) 
to the same scale. The dashed and dot-dot-dot-dashed lines 
show the degree of polarization for these two cases, normal- 
ized to the maximum polarization for that phase function, c. 
Model 3 (ellipsoidal cloud). The solid line shows the profile of 
the model (source plus scattered light), in which the scatter- 
ing cloud has an eccentricity e = 0.9, and the dot-dashed line 
shows the profile of Model 1 to the same scale. The dashed 
and dot-dot-dot-dashed lines show the degree of polarization 
for these two cases. 



Parameter 


Model 1 


Model 2 


Model 3 


^scat 


0.13 


0.13 


0.13 


9 


0.0 


0.6 


0.0 


e 


0.0 


0.0 


0.9 


Rc (") 


18 


44.5 


12.5, 28.5 


ijh ( cm 3 ) 


3300 


1200 


4000 


P 


0.35-0.7 


0.05-0.15 


0.55-0.75 



line). In each case, the source profile of a thin spherical 
shell (of radius 3") has been added (dotted line) and the 
scattering optical depth to the source has been taken to be 
0.13. The filled circles signify the measured intensity at 2" 
intervals along the slit and the horizontal scale has been 
adjusted to give the best fit between model and observa- 
tions. The degree of polarization is also shown. It can 
be seen that the scattered intensity hardly differs between 
the point and extended source, although the polarization 
is lower for the extended source. The cloud radius deduced 
from the fit is 72 c = 18 " , which, using r scat = 0.13, implies 
a number density for the upstream gas of «h — 3300 cm -3 . 

Model 2 is also for a spherical cloud but this time us- 
ing a modified H-G phase function with g = 0.6. Again, 
T scat is taken to be 0.13 and the Rayleigh scattering model 
with the same parameters is shown for comparison (dot- 
dashed line). The forward-throwing scattering causes the 
intensity to fall off more sharply with distance than in 
the Rayleigh case and so the horizontal scale needs to be 
stretched to agree with the observations, leading to a cloud 
radius i? c = 44.5" and a number density riH = 1200cm _3 . 

Model 3 reverts to a Rayleigh phase function but uses an 
ellipsoidal cloud (the major axis of which lies in the plane 
of the sky) of eccentricity e = 0.9 (a/b ~ 2.3). As usual, 
T scat = 0.13 and a spherical cloud with radius equal to the 
minor axis of the ellipsoid is plotted for comparison (dot- 
dashed line). In contrast to Model 2, the intensity now 
falls off less steeply with distance (and extends out beyond 
R c ) and so the horizontal scale needs to be compressed 
to fit the observations, leading to R c = 12.5", R c ' = 
28.5" and n H = 4000cm- 3 . 

The important features of the three models are summa- 
rized in Table 1. Also included in the table is a range of po- 
larizations predicted for the scattered light. The first value 
is for 4" upstream of the source and the second value is 
for 12" upstream. For non-Rayleigh scattering (Model 2), 
a Rayleigh form is still used for the polarization phase 



6 



function but with the maximimum degree of polarization 
Pmax (at a scattering angle of 90°) allowed to be less than 
unity. This is a reasonable approximation to experimental 
results for scattering from irregularly shaped particles (see 
Giese et al. 1978 and discussion in Paper III). The value 
of Pmax appropriate to the phase function used in Model 2 
is~0.5. 

Because it is the slope of the brightness profile that is 
being fitted here, the vertical positioning of the observa- 
tional points as a whole with respect to the model curve 
has been allowed to vary when finding the best fit. How- 
ever, for all three models it turns out that the ratio of 
scattered to source intensity in the slit is within 20% of 
that measured. Considering the uncertainties in this mea- 
surement, this is strong confirming evidence that the scat- 
tering optical depth adopted for the models, determined 
from the total scattered intensity, is essentially correct. As 
can be seen, the sizes, and hence densities, of the models 
vary considerably. Model 2 has the lowest density because 
its strongly forward-throwing scattering phase function re- 
quires a very large cloud to reproduce the observed bright- 
ness profile. Model 3, on the other hand, has to have a 
high density because its ellipsoidal shape means that it 
extends much further in the direction of the slit than per- 
pendicular to it and is hence smaller than a spheroidal 
cloud would be with a similar slit profile. Note that the 
predicted polarizations, also, have a wide spread between 
the models. Model 3 has the largest polarization because 
its elongated shape leads to the scatterers being concen- 
trated towards the plane of the sky, with scattering angles 
nearest 90°. Model 2 has the lowest polarization because, 
first, the high g means that most of the scattered light 
is scattered through small angles and, second, the peak 
polarization for this phase function is only 0.5. 

Since it is quite possible that asymmetric scattering and 
an aspherical scattering cloud are present in HH 1, then 
it could be argued that the extremes of Models 2 and 3 
will cancel each other out, leaving the parameters derived 
from Model 1, despite the naivete of its assumptions, as 
the most appropriate. However, since it is the dust with 
a projected position near the source that is most impor- 
tant for forward-throwing scattering, the eccentricity of 
the cloud is relatively unimportant for high values of g, as 
can be seen from Figure 2, and an ellipsoidal cloud would 
behave in a very similar fashion to a spherical cloud of ra- 
dius equal to its semiminor axis. Hence, Model 2 should be 
roughly appropriate to non-Rayleigh scattering, whatever 
the geometry of the cloud. 

Although Models 1, 2, and 3 all fit the observed dis- 
tribution of scattered light quite well, they also predict 
embarrassingly large densities for the scattering region, 
especially Models 1 and 3. Although line ratio diagnos- 



tics (Solf, Bohm & Raga 1988) imply electron densities of 
~ 3000-10000 cm -3 at the bow shock (where the gas will 
be fully ionized), this is gas that has just passed through 
a radiative shock and the pre-shock gas should be at a 
much lower density. The line ratios upstream of the bow 
shock are no guide to the density there, since most of the 
light is scattered, but Hartigan (1989) estimates nn~200- 
400 cm -3 , using general arguments from the kinematics 
and dynamics of the jet. This is much lower than that 
derived from the models, which suggests either that the 
material directly upstream of the bow shock is at a lower 
density than is typical of its vicinity or that too low a value 
of <7 H has been used. The latter option is quite possible 
since values of this quantity derived for Hll regions can 
vary by an order of magnitude (Osterbrock 1974). The 
first option, on the other hand, could pertain if HH 1 was 
moving into a channel that had been swept out by previous 
episodes of jet activity (see the discussion of time- varying 
jets in the introduction). 

2.3 Spectral Line Shapes of the Scattered 
Light 

The peak velocity and velocity width (FWHM) of the scat- 
tered [Sll] light, as a function of position along the slit, can 
be found from Figure 5 of SB. It is also possible to esti- 
mate the skewness of the line shapes by combining this 
information with that contained in the longslit spectro- 
gram. If the half- width-half-maxima to the blue and to 
the red of the peak velocity are denoted by Ai>bi U e an d 
Ai> re d respectively, then a simple estimate of the skewness 
is 

At; b iue-Ai; red 

K = ~^ a • (12) 

Ai; b iue + Ai; red 

These quantities (i> pe ak, A 

Dfwhmi and k) are shown m Ta- 
ble 2 for various positions along the slit. 

In Figure 4 these observational data are compared with 
predictions from the scattering models that successfully 
reproduced the spatial distribution of scattered light in 
the previous section. There are now two further relevant 
parameters for the models: the source velocity magnitude 
«o and the inclination of this velocity to the plane of the 
sky a. Since the intrinsic light from the bow shock of HH 1 
has an extremely low radial velocity, it is obvious that the 
motion of the HH 1 jet is very close to the plane of the sky. 
Hence, to start with, a = will be assumed. The following 
points should be borne in mind when interpreting these 
results: 

1. The observational velocities are all heliocentric, 
whereas the model velocities assume the scatterers 
have no velocity component towards the observer. 
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Figure 4: The velocity moments of the total scattered in- 
tensity (solid line) and the polarized intensity (dashed line) 
as a function of upstream distance from the source for the 
three models of Table 1. In each case, a = is assumed and 
the source velocity magnitude Mo is varied so as to provide 
the best fit to i> pe ak- The filled circles show the observa- 
tional results (total intensity) of SB. a. Model 1, Rayleigh 
scattering, R c = 18", u = 80 km s" 1 . b. Model 2, Mod- 
ified Henyey-Greenstein scattering, g = 0.6, R c = 44.5", 
«o = 270 km s _1 . c. Model 3, Rayleigh scattering, ellipsoidal 
dust cloud, e = 0.9, semimajor axis 28.5", Mo = 70kms _1 . 



TABLE 2 

Lineshape Parameters of the [Sii]AA6716, 6730A 
Emission Lines 



R — Rq 

(") 


^peak 

( km s _1 ) 


FWHM 

( km s _1 ) 


K 





0±9 


79±20 




2 


-36±9 


113±25 


0.46±0.12 


4 


-41±9 


150±25 


0.60±0.07 


6 


-74±6 


167±25 


0.24±0.10 


8 


-94±16 


148±40 


0.18±0.20 


10 


-119±19 


202±40 


0.04±0.12 



Note. — Measured at various distances upstream from the edge 
of the leading condensation in HH 1. All values were measured 
from the data of SB. The errors are estimated from the "nois- 
iness" of the data and do not take into account any contami- 
nation of the scattered light by intrinsic emission. 

This is probably unimportant since the line-of-sight 
velocity of the molecular cloud associated with HH 1 
is very low (Rodriguez 1989). 

2. The error bars on the observational points were esti- 
mated from the noisiness of the data at each position. 
They are much lower than the formal velocity resolu- 
tion of the spectrograph (55kms _1 ). 

3. The observational velocities were determined by fit- 
ting Gaussians to the line profiles, whereas, with the 
model velocities, the maximum of each line profile was 
determined. For a skew line profile, the difference be- 
tween the velocities deduced by these two methods 
could be of order ~ 0.5kAi> F whm- 

Also shown in the figures are the corresponding quantities 
for the polarized intensity. For each of Models 1,2, and 3, 
«o is varied so as to find the best fit between the observa- 
tions and model predictions for i> pe ak- It can be seen that, 
for Models 1 and 3, the fit is so poor that the choice of Mo 
is quite arbitrary. The poor fit is because the observations 
show a gradual increase in the blue shift of the light as 
the distance from the source is increased, whereas Mod- 
els 1 and 3 show a roughly constant blue shift of — Mo for all 
upstream positions (except for very close to the source in 
Model 1). This can also be seen in the relevant position- 
velocity diagram (left panel of Figure 5), where the scat- 
tered intensity shows a ridge along v = — 1 for x greater 
than the source radius (zero in the case of Model 3). It 
is also clear from the further position- velocity diagrams in 
Figure 13 of Paper I that varying a will not improve the 
quality of the fits, since the ridge still exists for a = 0.5. 
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a = 0.0, Aperture A, g= 0.0 

Lowest contour = 3.5 x 10" 4 




a = 0.0, Aperture A, g= 0.6 

Lowest contour = 2.4 x 10" 4 
Maximum polariza tion I / - 

1 " 



near source 



far from source 




■1 -0.5 0.5 
Velocity Shift 



-1 -0.5 0.5 1 
Velocity Shift 



Figure 5: Position-velocity diagrams of scattered light cor- 
responding to two different models for the dust scattering 
phase function. Left Panel (from Figure 13 of Paper I): 
Rayleigh scattering (Model 1). Note that the peak inten- 
sity of the scattered light is at v = —l (corresponding to the 
source velocity) for all upstream positions. Right Panel: 
Non-Rayleigh scattering (Model 2). Note that the peak in- 
tensity of the scattered light is increasingly blue-shifted with 
increasing distance upstream of the source. 




v r km s" 1 v r km s" 1 

Figure 6: Schematic illustration of how contamination by 
intrinsic emission can affect the measured line widths. The 
scattered component is shown by the dashed line and the 
intrinsic component by the dotted line. The width of the 
total profile (solid line) is smaller near the source than it 
is far away from it, despite the fact that the width of the 
scattered component does not change. 



Model 2, on the other hand, fits the observed velocities 
much better. This is because the asymmetry of the scat- 
tering causes the majority of the scattered light to come 
from dust on the near side of the cloud, which sees an in- 
creasingly blue-shifted source as the distance upstream is 
increased, whereas, with Rayleigh scattering, it is the dust 
directly upstream of the source which contributes most 
to the scattered intensity from a given line of sight and 
this dust sees the full blue shift. The phenomenon can 
also clearly be seen in the position-velocity diagram (right 
panel of Figure 5), where a ridge in the scattered intensity 
roughly follows the line v = —x. Even Model 2, though, is 
not a perfect fit to the observations of i> pe ak, particularly 
near the source, where the observations suggest i> pea k —> 0. 
This, however, may be due to contamination by intrinsic 
emission. 

If the velocity width Ai> F whm is now considered, it can 
be seen that, again, Models 1 and 3 are very bad fits 
to the observations. They both show Ai> F whm decreas- 
ing sharply with distance, whereas the observations show 
it increasing. Model 2 shows a roughly constant Ai> F whm, 
which also appears to be inconsistent with the observa- 
tions. However, as illustrated in Figure 6, contamination 
by intrinsic emission would lead to A 

Dfwhm being under- 
estimated near the source and yet overestimated far from 
the source. This apparent paradox arises in the following 
manner. The scattered and intrinsic light in the vicinity of 
the source have similar velocities but the scattered light is 
much broader so that, even if the intrinsic light is fainter 
than the scattered light, the FWHM of the sum of the two 



components is essentially that of the intrinsic light. Fur- 
ther upstream, however, the intrinsic light, which is prob- 
ably from a foreground nebula and which can be clearly 
seen in the Ha spectrogram of SB, is much less blue shifted 
than the scattered light so that the FWHM of the blend of 
the two components is larger than that of either of them 
individually. If this hypothesis were true, then it may 
also affect the peak velocity of the light at large distances. 
There is slight evidence that i> pea k starts to increase again 
beyond ~ 12" from the source but the data here are very 
noisy. There are also faint narrow knots of emission visi- 
ble near zero velocity from 16 to 26" upstream in the [Sll] 
spectrogram. 

Turning now to the skewness k of the line profiles, it can 
be seen that, yet again, Model 2 agrees, at least qualita- 
tively, with the observational results, while Models 1 and 3 
do not even produce the correct sign. Since the skewness 
is, in effect, the third velocity moment of the specific in- 
tensity and was very crudely determined from the data, it 
is heartening that the agreement with Model 2 is as good 
as it is, especially since k will also be very sensitive to in- 
trinsic contamination (which perhaps explains the decline 
at large distances from the source). 

The effects of extinction on the line profile of the total 
scattered light from the cloud are similar in some ways to 
those of asymmetric scattering, as will be demonstrated 
in Paper III, so it might be thought that the inclusion 
of extinction would allow Model 1 to fit the observations 
better. This is not the case, however, since even using an 
extinction optical depth r ext = 1.0 does not significantly 
change i> pea k of the upstream scattered light. This can 
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Figure 7: The same as Figure 4 but for Model la, which dif- 
fers from Model 1 in that extinction has been included with 
an optical depth r ext = 1.0. The source velocity magnitude 
uo = 90 kms -1 and a = 0. The filled circles show the ob- 
servational results (total intensity) of SB and it can be seen 
that even such a large optical depth does not improve the fit 
significantly. 
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Figure 8: The same as Figure 4 but for Models 2c and 2d, 
which differ from Model 2 in that they have different in- 
clinations of the source velocity to the plane of the sky a. 
a. Model 2c, a = -O.lrad ~ -6°, u = 330kms -1 . b. 
Model 2d, a = -0.2rad~-12°, u = 560 km s -1 . The filled 
circles show the observational results (total intensity) of SB. 



TABLE 3 

Source Velocity m from Scattering Models 



Model 


9 


e 


Text 


a 


MO 


RMS dev 










(°) 




[ km s -1 ) 


1 


0.0 


0.0 


0.0 


0.0 


80 


28 


la 


0.0 


0.0 


1.0 


0.0 


90 


26 


2 


0.6 


0.0 


0.0 


0.0 


270 


15 


2a 


0.6 


0.0 


0.0 


5.7 


190 


22 


2b 


0.6 


0.0 


0.0 


11.5 


150 


27 


2c 


0.6 


0.0 


0.0 


-5.7 


330 


14 


2d 


0.6 


0.0 


0.0 


-11.5 


560 


10 


3 


0.0 


0.9 


0.0 


0.0 


70 


38 



Note. — Derived from an unweighted least-squares fit to the 
observed values of the blue shift of the scattered light « pe ak- 
Also shown is the RMS deviation of the observational data 
from this fit. 



be seen in Figure 7, in which the velocity moments are 
plotted for such a model (Model la). 

Accepting that g = 0.Q provides the best fit to the obser- 
vations, the effects of varying a will now be investigated, 
to see if the fit can be improved. Table 3 summarizes the 
parameters and quality of fit for all the models considered 
so far, together with four further models (2a-2d), which 



are identical to Model 2 except for different values of a. It 
can be seen from the table that Models 2a and 2b (positive 
a) fit the data less well than Model 2. Models 2c and 2d 
(negative a), on the other hand, seem, at first glance, to 
fit the data better and are shown in Figure 8. However, 
the source velocity deduced from the fit of Model 2d is 
560 kms -1 , which is far too high. Raga, Barnes & Ma- 
teo (1990) found a velocity for HH1 of 350 kms -1 from 
proper motion measurements and it seems highly unlikely 
that «o can be higher than this. This model would also 
predict that the direct light from the source should be red- 
shifted by «o sin 0.2 ~ 100 kms -1 , which is not observed, 
and produces values of Ai> F whm much higher than those 
seen. Model 2c predicts Mo = 330 km s -1 , which is possible 
if the upstream dust is stationary, but this model is also 
rather inconsistent with the observed direct light from the 
source, which suggests a bow shock moving in the plane of 
the sky. Model 2c does, however, produce the best fit to 
the observed A 

Dfwhmi although this is quite sensitive to 
the assumed velocity width of the source. If a higher value 
had been used (perhaps more realistically) then Model 2 
would fit the observed A 

Dpwhm best. 
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3 Discussion 

In conclusion, Models 2 and 2c have been most successful 
in reproducing the velocity characteristics of the scattered 
upstream light in HH 1. They differ most significantly in 
their implied values for the magnitude of the relative ve- 
locity between the source and the scatterers uq. Model 2 
predicts uq = 270 km s _1 , which, although still somewhat 
higher than the ~ 200kms _1 found by fitting bow shock 
models to the line ratios and line profiles of HH 1 (Raga 
et al. 1988), is lower than that derived from proper mo- 
tion measurements. Hence, if the proper motion values 
are taken at face value to represent the true speed of the 
shock-excited gas in HH 1 (rather than merely a pattern 
speed), then Model 2 provides some support for the no- 
tion that the bow shock of HH 1 is propagating into a 
moving environment. Such a situation could arise as the 
result of a variable velocity jet, as discussed in the intro- 
duction. Model 2c, however, lends no support to such a 
notion since it predicts a uo similar to the proper motion. 
Although current observational evidence makes it difficult 
to discriminate between these two models, spectropolari- 
metric observations of the scattered light could help con- 
siderably. This is because, as can be seen from the figures, 
the peak velocity of the polarized intensity is roughly con- 
stant at i> P eak ~ 0.9«o for all upstream displacements and 
for all the models. Measurement of this quantity would 
hence provide a relatively unambiguous determination of 
«o- The one factor that might spoil this rosy picture is 
the uncertainty over the source line profile (as seen by the 
scatterers). All the above models were calculated using a 
Gaussian profile for the source with a dimensionless veloc- 
ity width w = . 1 , corresponding to a F WHM of 45 km s~ 1 
for Model 2, or 55kms _1 for Model 2c. This is slightly 
lower than the observed value for the source (79kms _1 ) 
but theoretical models of bow shock spectra (Raga 1985) 
show that, although, for a bow shock in the plane of the 
sky, many emission lines have velocity widths compara- 
ble with the velocity of the bow shock, the upstream dust 
would see a much narrower line profile, albeit with broad 
faint wings. Figure 9 shows the result of Model 2e, which 
is identical to Model 2 but with id = 0.2. For this model, 
the best-fit value of uo is 230kms _1 , which is much more 
in keeping with the bow shock models and the moving 
medium hypothesis. The RMS deviation of the fit from the 
observed velocity peaks is, at 17kms _1 , only marginally 
worse than Model 2 and the fit to Ai> F whm is considerably 
better. The polarized intensity i> pea k is at~0.85«o for this 
model, so it can be seen that if the source velocity width 
is uncertain by a factor two, there will be an uncertainty 
in the derivation of uo by this method of ^ 10%. 

Although the results presented here seem to rule out a 
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Figure 9: The same as Figure 4 but for Model 2e, which 
differs from Model 2 in having a source line profile with dou- 
ble the dimensionless velocity width. The source velocity 
magnitude «o = 230 km s -1 . The filled circles show the ob- 
servational results (total intensity) of SB. 



Rayleigh form for the scattering phase function, it is quite 
possible that such a model may be capable of being fitted 
to the observations by, for example, adopting a density 
distribution different from the homogeneous one assumed 
here or by introducing an internal velocity field among 
the scatterers. However, such further complications would 
necessarily be rather ad hoc although their introduction 
may become necessary, once spectropolarimetric observa- 
tions are obtained, if the simple models developed so far 
no longer suffice. 

To summarize, the results of this case study demon- 
strate the powerful potential of emission line spectropo- 
larimetry for studying astrophysical flows. The relatively 
simple geometry of HH 1 allows the polarized scattered 
light to reveal an "image" of the bowshock at the head 
of the jet, as seen from the upstream dust. The Doppler 
shift of this light is then a direct measure of the velocity 
of the bowshock with respect to the dust. It is therefore 
of paramount importance to conduct spectropolarimetric 
observations of this object in order to constrain the vari- 
ous possible models. 
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